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1. Introduction 

An accurate determination of the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements, Vij, 
is an essential step in testing the Standard Model (SM). Like the quark masses, the couplings 
Vij are free parameters of the SM and therefore cannot be predicted. Instead they are extracted 
after confronting the experimental measurements to the SM theoretical expressions. The simplest 
processes in that respect are the leptonic and semileptonic decays of pseudoscalar mesons. In this 
note we consider D decays, namely 

( D ^ V v ( ) = |vy2_^A-VV)l^)l 2 , 

T(D+ -> tV t ) = \V cq \ 2 ^f D m D ml f 1 - j , (1.1) 

where I is \i or e. The left-hand-side in the above expressions is measured experimentally, while the 
computation of hadronic form factor, F + (q 2 ), and/or the meson decay constant, fp , requires a first 
principle description of non-perturbative QCD effects. We now restrain our attention to F + (q 2 ), one 
of the two form factors which parameterise the SM weak matrix element (n(k)\(V —A)^\D(p)} = 
(n(k)\V^\D(p)),i.e., 

(K(k)\V,\D(p)) = (p + k-q^^) F + (q 2 )+q,^^F (q 2 ), (1.2) 

both depending on q 2 = (p — k) 2 only, with q 2 £ (0, (mo — mx) 2 ]. 

Lattice QCD is the only currently available method which allows us to compute this matrix 
element without introducing any extra parameter and, at least in principle, with an accuracy that 
can be matched to the experimental one. In practice, however, there is still quite a room for im- 
provement on systematic errors. Here we want to address those that arise from the extraction of the 



matrix element ( |1.2| ) from the correlation functions. 

1.1 An abridged description of the standard procedure 

The standard method consists in computing the 2- and 3-point functions, namely, 

X 71 



cf" D (k,q;t,t SO urce) = £ {{qy5q)xu ur cA^ c )iAcy^)x),Q) e 

Q«*<w ra ^e-^(n(%)\V^\Dip))^e- E ^ "0 , (1.3) 

where we also indicate their asymptotic behavior, using the standard notation, = \ {0\cY^q\D{p)) | 2 , 
and similar for 3£ % . The matrix element ( |L2| ) corresponds to a plateau of the ratio 

R = __^3 ^ Source) ^ y °« f -^»- (n(k)\V^\D(p)) . (1.4) 

C2 {k^t source 



-i{qz-kx) 
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Figure 1: The valence quark diagram of the 3-point function, (k,q;t,t source ). We indicate the twisting 

angles, 61,2,3, which are discussed in the text. 



In case of D-decays the plateaus are known not to be long enough to guarantee a percent accu- 
racy (especially when the momenta are given to either of the two mesons). Furthermore, V ' 3? n 
and \J 2?d should be computed from a separate study of the 2-point functions, the errors of which 
are carried over to the ratio R. Finally, a non-negligible statistical error is introduced by the mul- 
tiplicative renormalization of the vector current, especially when consistently implementing the 
0(a) improvement of the Wilson quark operators on the realistic lattices and one of the quark be- 
ing charmed (heavy). All these difficulties can be avoided by considering various double ratios of 
3-point correlation functions. The efficiency of the double ratios was introduced and tested first 
in heavy-to-heavy ^ and then in light-to-light decays [g]. In D-decays to a light meson an extra 
problem is related to the available kinematics on the lattice. More specifically, with the periodic 
boundary conditions the minimal momentum on realistic lattices is (2k) /L, which is too large 
(L = Nia not large enough) if one is to keep the lattice spacing, a, sufficiently small in order to 
accommodate the charm quark mass. To get around this problem, we adopted the twisted boundary 
conditions (twBC) recently proposed in ref. [J3j] . In such a way the momenta p and/or k become 
6/L = (do, 60, 60) /L, where the components 60 are chosen anywhere between < % < 71. In this 
way we are able to compute the form factor F + (q 2 ) at several q 2 > 0. 

2. Double ratios 

In what follows we consider 4 different strategies to increase the accuracy of extraction the 
form factor F + (q 2 ) for several values of q 2 by combining the 3-point functions in suitable double 
ratios. At this point we should emphasize that in each of the strategies discussed in this section, 
the multiplicative renormalization factors, as well as the source terms, 2£ n jy, cancel out. We then 
numerically test each of the proposed strategies to check whether or not a plateau region is pro- 
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nounced enough and the statistical quality of the signal satisfactory to reach a percent accuracy of 
the extracted form factor. 

2.1 First strategy 

We first keep the D-meson at rest and inject momenta to the pion only. This is done by 
imposing d\ = 63 = (c.f. fig. [I]), while for 6 2 we choose several different values to explore the 
kinematics available from this decay, < q 2 < q 2 max . Similar to what has been proposed by JLQCD 
in their study of the -decay [|J], we consider the following double ratios (k = 2 /L): 
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These ratios then can be cast into the expressions leading to F + (q 

mo + m n 



1, i.e., 



F+tf) =R l xF Q (qi 
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(2.4) 
(2.5) 
(2.6) 



For short we wrote q 2 ^ (q 2 ) = (mj) — m 2 z )[Fo(q 2 ) /F + (q 2 ) — 1]. The quality of the plateaus is pre- 
sented in fig. ^. To that end we use the publicly available ensembles of the SU(3) gauge field 
configurations, produced by the QCDSF collaboration by using the & (a) -improved Wilson quark 
action with Np = 2 [|0. We computed the quark propagators and correlation functions on the con- 
figurations gathered at j8 = 5.29, corresponding to a ~ 0.08 fm, on the 24 3 x 48 lattice, each time 
keeping the light valence quark mass equal to that of the sea quark. In fig. ^ we see that the signals 
for Rq and R\ are indeed very good. The fact that the signal for R 2 is not as good does not trouble 
the whole strategy because it is only needed to compute ^(q 2 ), which itself is a correction to 1 in 



both eq. (16) and in eq. (2A). On the other hand the fact that the quality is less good for R 2 (t) is 
expected as its computation involves the correlation function with with "^"-matrices (mixing the 
"large" and "small" components of the Dirac spinors), in contrast to /?o,i(0 m which we compute 
the correlation functions with "70 "-matrix only. 

2.2 Second and third strategies 

Next we consider the kinematics in which we either keep the D-meson at rest (8\ = 63 = 0) 
and inject momenta to the pion source (d 2 = d q ), or keep the pion at rest (§1 = 63 = 0) and inject 
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Figure 2: Double ratios, Z?o,i,2(0> illustrating the first strategy. The lattice data refer to the simulation with 
Wilson G[a) -improved quarks with Np = 2 (j3 = 5.29) and Kg = K" sea = 0.1355. The shown signals refer to 
q 2 « 1 GeV 2 . 



momenta to D-meson (d\ = 9 C ). Notice that keeping the q 2 fixed requires that B c = (mj)/m K ) x 6 q . 
We then build the following four double ratios: 
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where p = B\/L, and k = B2/L. In terms of form factors, the above ratios read 

[F + (q 2 )] 2 , Ra 



[E D + m }t + (E D -m 7t )^(q 2 )] 2 2 2 _ [m D +E n + (m D -E n )% (q 2 )] 2 



2\i2 



2m K 2E D 



PL 



2E Jl 2m D 
^-^iF^ 2 )?. 



(2.8) 



Am n 1 ~ rK1 " •* • 4 mD 

Note that UKQCD ^ recently considered R3 and R' 3 to probe Ff^^fi). We now need to combine 
two of the above four double ratios to obtain F + (q 2 ), i.e., either R3 with R' 3 , R 4 with R' 4 , R3 with 
R' 4 , or R' 3 with /?4. In fig. |5| we show the quality of the signals corresponding to the same set-up 
as the signals displayed in fig. [2[ We observe that only R 4 and R' 4 are reasonably good, whereas 
/?3(?) and R' 3 (t) do not exhibit plateaus, likely due to the fact that mvlm n is large, so that \Q C \ in 
@c = {wiD/inn) x Qq is such that the twBC simply destroys the signal. Therefore, our tests suggest 
that the ratios R 4 and R' 4 can be used to compute the form factor F + {q 2 ) for various values of q 2 . 
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Figure 3: Illustration of the signals for the double ratios used in the second, third (R3, R4, R ! 3 , R' 4 ) and fourth 
strategies {R$, R(,), defined in the text. 



2.3 Fourth strategy 

If one wants to study the shape of the form factor (i.e., its q 2 -dependence), then it is worth 
trying to impose the twBC on the spectator (63 / 0) without twisting the other two quarks (di = 
62 = 0). In such a way we may probe many values of q 2 but never reach q 2 = 0, because in this 



set-up the condition Ed = E n does not allow a real solution in \p\ 
we consider the following two ratios: 



\k\ = 63 /L. To test this option 
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where index "0" is used to distinguish that the spectator quark is actually twisted. Expressed in 
terms of form factors 
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The corresponding signals from our numerical study are shown in fig. |3[ While the statistical 
quality of Rs(t) is reasonably good, the signal for Re(t) is not promising if we are after a strategy 
that could lead us to a percent accuracy on the extracted form factor value. Although we did not 
try it, we suspect the flatness of R5 could be achieved by a judicious choice of smearing. We point 
out, however, that our numerical evaluation of the ratios ( j2~^ ) indicate the large statistical errors so 
that this strategy is not competitive with the first or the third ones discussed in this section. 



3. Summary 

It this note we report on the results of our exploratory study in which we use various double 
ratios and twisted boundary condition on the quark propagators in order to extract the form factor 
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D -» Plv Form Factors - k - 0.1355 D -> Plv Form Factors - ic = 0.1359 




q 2 in GeV 2 q 2 in GeV 2 

Figure 4: D — > n semileptonic form factors for the unphysically heavy pions: left panel corresponds to 
m % ~ 750MeV, the right one to m n ~ 600MeV 



relevant to the semileptonic heavy-to-light D-decays to a percent accuracy. In total we proposed 
four different strategies, which we then tested numerically on the set of unquenched (Np = 2) gauge 
field configurations. On the basis of our analysis we conclude that the first (double ratios /?o,l,2) an d 
the third strategy (R4 and R' 4 ), discussed in the text, can be used to compute the form factor F + {q 2 ) 
to a desired precision. Even though our numerical tests are made by using the Wilson quarks, our 
conclusions apply to any lattice QCD action. 

In fig. |], we illustrate the results obtained by employing the first strategy at three values of the 
twisting angle 6. Those will be improved and the results discussed in our forthcoming paper. 
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